One-dimensional fermionic gases with attractive p-wave interaction in a hard-wall trap 
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We investigate the ground state of the one-dimensional fermionic system enclosed in a hard- wall 
trap with attractive contact p-wave interactions. Based on the Bethe ansatz method, the explicit 
wave function is derived by numerically solving the Bethe ansatz equations for the full physical 
regimes (— oo < cf < 0). With the exact wave function some quantities which are important 
in many-body physics are obtained, including the one-body density matrix and the momentum 
distribution of the ground state for finite system. It is shown that the shell structure of the density 
profiles disappears with the increase of the interaction and in the fermionic Tonks-Girardeau (FTG) 
limit the density distribution shows the same behavior as that of an ideal Bose gas. However the 
one-body density matrix and the momentum distribution exhibit completely different structures 
compared with their bosonic counterparts. 
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I. INTRODUCTION 

Experimental realization of trapped one-dimensional 
(ID) cold atom systems [l|, [H, H E| are triggering more 
and more theoretical efforts to study the ID many-body 
physics beyond the mean-field theory. For the ultracold 
quantum gases tightly confined in waveguides, the dy- 
namics are effectively described by a ID model due to the 
radial degrees of freedom are frozen (HQ. Further, the 
ability of tuning the effective ID interactions by Feshbach 
resonance allows experimental access to the very strongly 
interactingregime where correlation effects are greatly 
enhanced 0, [1, 1^. [lfjl. Hlj . In the limit of the Tonks- 
Girardeau (TG) PJl] gas with effective coupling constant 
giD — > oo, the many-body problem of a TG gas can be 
mapped to that of a free Fermi gas by the Bose-Fermi 
mapping, which has been verified by two experimental 
groups [H, H[ . This Bose-Fermi duality was generalized to 
show the equivalence between a ID fermionic system and 
a bosonic one with the reversed role of strong and weak 
couplings [13J. Recently, the exact ground state of the 
fermionic TG (FTG) gas, defined as a ID spin-polarized 
fermionic gas with infinitely strong attractive p-wave in- 
teractions, has been determined by inversely Fermi-Bose 
mapping to the ideal Bose gas Q [H, [H, [TtJ • 

A key experimental challenge is to obtain superfluidity 
with pairs in nonzero orbital angular momentum states 
by using p-wave, or maybe even d-wave Feshbach reso- 
nances. In general, the p-wave interaction is very weak 
comparing with the s -wave interaction. However, for 
a spin-polarized fermionic gas, the s-wave scattering is 
forbidden due to the Pauli exclusion principle and thus 
the p -wave interaction is dominant. Furthermore, the 



p-wave interaction can be g reatly enhanced by the Fesh- 
bach resonances [H, [H, H3, H3 and using a p-wave Fesh- 
bach resonance between 40 K atoms Jin's group at JILA 
have successfully produced and detected molecules with 
lifetimes on the order of milliseconds on both the BEC 
and the BCS side of the resonance [22j. For a ID gas, the 
additional confinement induced resonance permits one to 
tune the ID effective interaction via a 3D Feshbach res- 
onance @, H, HU . 
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In this paper, we report on a detailed study of the ID 
Fermi gases in the infinitely deep square potential well. 
We will show that the model of fermionic gases with at- 
tractive p-wave interactions in such a one-dimensional 
hard-wall trap is exactly solvable by the Bethe-ansatz 
method. The experimental efforts in trapping ultra- 
cold gases near micro-fabricated surfaces, the so-called 
"atom chips" [23|, [24j], and various innovative features 
in designing the optical box trap [2^], are specifically 
aimed at studying the surprisingly rich variety of physical 
regimes predicted for the ID Bose gas and have stimu- 
lated many theoretical studies on the physics in a box 
trap [26|, [23, l28j |- Different from the harmonic trap, the 
interacting model in a hard-wall trap is integrable and 
thus could provide us some exact pictures for understand- 
ing the trapped many-body systems. So far, there has 
been a growing interest in the exactly solved models in 
the hard- wall trap [26|, [2^, l3(ij |. but most of them fo- 
cus on the Bose gas and the Fermi gas with odd-wave 
interactions is not addressed. While the theoretical un- 
derstanding of the correlation effect of bosonic system 
has been investigated extensively [EBBS EH El > the 
fermionic system is not well understood except in the 
so-called FTG limit EE El El . 
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II. FORMULATION OF THE MODEL AND ITS 
EXACT SOLUTION 

We consider an TV-particle system with finite, attrac- 
tive p-wave interaction in a one-dimensional box of length 
L, which obviously fills the gap between the FTG limit 
and free Fermions. The Schrodinger equation can be for- 
mulated as 



N 

£ 

i=l 



K_cP_ 

2m dx 2 



l<i<j<N 



* = (1) 



where V (xi — Xj) is the pseudo-potential describing the p 
-wave scattering. It has been shown that the p-wave scat- 
tering of two spin-polarized fermions in a tightly confined 
wa veg uide can be well described by the contact condition 



*f (a 



-af D — ^ F (x l = x J ±0),(2) 



where 



1 1D 



3C(3/2) fop 
2V2tt \l± 



(3) 



is the effective ID scattering length with a p the p-wave 
scattering length and l± — y^h/mu)± the transverse os- 
cillator length [l4j| . The contact condition can be repro- 
duced by using the following pseudo-potential [3, [32] 



V(x) 



m ox ox 



(4) 



where x = xi — Xj and d x — (d Xi —d Xj )/2. The scattering 
length can be tuned readily from to — oo by sweeping 
an external magnetic field - the Feshbach resonance, or 
by changing the geometry of the trapping potential - the 
confinement induced resonance and in this paper the full 
physical regimes — oo < af D < will be studied. Similar 
to the case of Bose gas, the important parameter char- 
acterizing the different physical regimes of the ID Fermi 
gas is 7 = mgf D p/h 2 , where gi D — —2Ti , a\ D jm and 
p = N/L. 

A standard rescaling procedure brings the Schrodinger 
equation into a dimensionless one (for simplicity we keep 
the original notations) 



with 
H = 



HV (x u --- ,x N )=EV(xx,--- ,x N ) 



v dxi dxj 



i—l 1 i<j 

where in the dimensionless interaction constant 2cf = 
\J 2m/h 2 a^ D we intentionally keep the factor 2 in accor- 
dance with the bosonic case 12 a . The wave function takes 



the general form 

* (Xl, • • • ,X N ) = ^0 (x qN - X qN _ 1 

Q 

X IfQ (Xq 1 , Xq 2 , ' ' ' , Xq N ) . 



j • • • 6 (Xq 2 X qi ) 



(5) 



where we have used Q to label the region < x qi < x q2 < 
• • • < x qN < L. The wave function of Fermions should 
follow the antisymmetry of exchange, so our model is 
simplified into the solution of 

H <PQ {%qi ,Xq N ) = EtfQ ( Xqi , • • • , X qN ) (6) 

in the region Q with the open boundary condition 

¥>Q (0, Xg 2 , ■ ■ ■ ,X qN )=(fiQ (x qi , x q2 , • ■ ■ , L) = 0. 

Using the Bethe ansatz method we obtain the wave- 
function parameterized by the set of quantum number 
k\ , k% , • • • , /cat known as quasi-momenta or rapidities [2r3 | 



'PQ i^qi ' x q2 ' ' ' ' j ZqN ) 
= (-l) Q ^(-l) P ^ P exp 



'iV-l 
Kj 



sin (kp 1 Xq 1 ) 



>< n 

Kj<JV 



Kj 

x exp (ik PN L) sin (k PN (L - 



with uj ab = arctan £[% b ~ arctan k b +k a and A PlP2 ... PN = 

llj<i (ikpi ~ ikpj ~ 1 / c f) {ik P i + ^Pj ~ V c f)- Here 
(— 1)** = ±1 and (— l) p = ±1 denote sign factors associ- 
ated with even/odd permutations of Q = (<j"i, q2, ■ ■ ■ , ?jv) 
and P = [px,Pi, ■ ■ • ,Pn)i respectively. The quasi- 
momenta fei, , • • • , fcjv can be easily determined from 
the Bethe ansatz equations 

t r kj - ki - i/c F kj + h - i/c F 
exp(i2kjL)= || ; 



i=H¥=j) 



ki + i/cp kj + ki + i/cp 



with j = 1, 2, • • • , N. These quasi-momenta lead us im- 
mediately to important physical quantities for our sys- 
tem. For example, the energy eigenvalue is given by 



N 



E = J2 k 



(7) 



and the total momentum by K — Y]j—i kj. It is clear 
that these Bethe ansatz equations are the same as those 
in the case of Bose gas if we simply make a substitu- 
tion c = —1/cp [26|. In the regime Q, there is one-to- 
one correspondence between the quasi-momentum solu- 
tion of attractive p-wave Fermi gas and that of repulsive 
Bose gas, but the total wave functions ^f(xi, X2, ■■ ■ ,Xn) 
take different forms because of their distinct exchange 
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symmetries. This difference can be easily seen from the 
one-body density matrix and the momentum distribu- 
tion. Taking the logarithm of Bethe ansatz equations, 
we have 

/ 1/cf 1/cf \ 
K,L = n,7r + > arctan arctan 

(8) 

For the ground state the set of integer rij = 1 (1 < 
j < N). For simplicity, in the following evaluation the 
length L will be taken as one. It is obvious that the solu- 
tions of kj is only relevant to cf/L for different L. The 
subsequent procedure is standard. By numerically solv- 
ing the transcendental equations eqs.©, we obtain the 
quasi-momenta kj and thus the ground state wave func- 
tion. In principle, all necessary information about the 
system can be inferred, including the one-body density 
matrix, the momentum distribution, and the excitation 
spectrum. Furthermore, one can apply the thermody- 
namic formalism for dealing with the one dimensional 
interacting systems developed by Yang and Yang [34j . 

Before we proceed to the general case with interme- 
diate interaction strength, we'd like to take a look at 
the situation in the two limiting cases. In the limit of 
free Fermion case, we have kj — jir (j = 1, • • • , N). In 
this exceptional case of cf — 0, the quasi-momentum ac- 
quires the genuine physical meaning of momentum and 
the atoms occupy the single particle momentum states 
according to Pauli exclusion principle. We thus have the 



total wave function 

* (xi,X2, ■ ■ ■ ,xn) 

Q 

N 

P 3=1 

with C the normalization constant, which is nothing 
but the Slater determinant of the lowest TV eigenstates 
sm(jirx) of the system of single particle in the hard 
wall potential. In the other limit of strongly attrac- 
tive interaction, i.e., the FTG limit, all quasi-momenta 
take the same value kj — ir (j — 1, • • • ,N) and 

(p Q (x qi ,x q2> --- ,x qN ) = C(-l) Q n^Li sin7r: %- Tllus 
we have the total wave function 

* (xi,x 2 , • • • ,x N ) 

N 

= C^2(-l) Q 8(xq N -z<w-i) ■■■S(x q2 -x qi ) JJsin7rx 9r 

Q 3=1 

On the other hand, the Fermi-Bose mapping method has 
been used to give exactly the ground state of the FTG 
gas in 

N 

* (x±,X2, ■ ■ ■ ,xn) = CA (xi,X2, ■ ■ ■ ,Xn) Y\_ siniTXj, 

3 = 1 
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with the antisymmetric function A (x\,X2, • •• , %n) = 
rii<j</<Ar s g rL ( a; i — x j)- ^ ls obvious that the above two 
wave functions match each other. 



III. GROUND STATE PROPERTIES AND 
COMPARISON WITH BOSONS 

We now turn to the system with finite interaction 
strength — oo < cp < 0. In this case the quasi-momenta 
are decided by numerically solving the Bethe ansatz 
equations and the total wave function is obtained by 
Eq. (O through ifQ {x qi , x q2 , • • • , x qN ) under the restric- 
tion of exchange antisymmetry. For the one dimensional 
interacting system, a quantity of fundamental impor- 
tance in many-body physics is the one-body density ma- 
trix, which, in terms of the ground state wave function 
\l> (xi, ■ ■ ■ ,Xff), is given by 

p(x,x') 

_ N / Q L dx 2 ■ ■ ■ dxN^* {x,X2, ■ • • , Xn) * {x',X2, ■■■ , x N ) 
J L dxi ■ ■ -dx N |* (xi,x 2 ,- ■ ■ ,x N )\ 2 

This quantity furnishes the expectation values of sin- 
gle particle observables such as the position density dis- 
tribution p(x) = p(x, x')\ x =x' j and the momentum dis- 
tribution which is simply the Fourier transformation of 
p(x,x'), 

n(k) = — ( L dx [ L dx'p(x, x')e- ik i x - x ') . (9) 
2?r Jo Jo 

We display the one-body density matrix and the po- 
sition density distribution in Fig. 1 and in Fig. 2 for 
different p-wave attractive interactions. The one-body 
density matrix expresses the self correlation and it means 
the probability that two successive measurements, one 
immediately following the other, will find the particle at 
the point x and x' , respectively. We notice that for all in- 
teracting strengthes there exists a strong enhancement of 
the diagonal contribution p(x, x') along the line x = x' . 
To see this more clearly the position density distribu- 
tion p(x) for N — 4 particles is shown in Fig. 2 with 
more variable interacting strengthes. For noninteracting 
system the atoms behave as ideal Fermions and the den- 
sity profiles show obvious spatial oscillation structure. 
Increasing the strength of the attractive p-wave interac- 
tion, as seen in Fig. 2, first leads to the depression of 
the amplitude of the oscillations of the density profile, 
followed by the emergence of the typical Gaussian-like 
bosonic behavior and contraction of the half- width of the 
density. In the limit of infinitely strong attractive inter- 
action between the Fermions, the system enters into the 
FTG regime, and the density shows the same smooth 
profile as that of noninteracting Bosons. Particularly, 
the density distribution in Fig. 2 for FTG is identical 
to that of the free Bose gas and there exists one-to-one 
correspondence for the density distribution between the 
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FIG. 3: (color online) Momentum distribution of Fermions 
with p-wave attractive interactions for N — 4. 
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FIG. 4: (color online) Momentum distribution of Bosons with 
repulsive interactions for N = 4. 



attractive fermionic gas and the repulsive bosonic gas 
with the interacting strength related by cf — — 1/c [261 ]. 

It is worth to emphasize that the density distribution is 
identical for TG bosons and noninteracting fermions, and 
also for FTG fermions and noninteracting bosons. Never- 
theless, their momentum distributions show remarkable 
differences as shown in Fig. 3 for fermions with different 
attractive interaction strengths and in Fig. 4 for bosons 
with the corresponding repulsive interaction strengths. 
The momentum distribution for fermions oscillate in the 
full regime and the number of oscillation peaks remains 
equal to the number of atoms in the system. It becomes 
more and more nonuniform in the momentum space and 
and two sharp spikes indicating that there is high prob- 
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FIG. 5: Gray scale plot of the one-body density matrix of 
Bosons with repulsive interactions for N = 4. (a) c = 0; (b) 
c = 10; (c) c = 50; (d) c = 100; (e) c = 1000; (f) c = oo. 

ability of finding the atom in momentum states around 
k *~ 1. Two other peaks at higher momentum diminish 
with increasing interaction but remain there even in the 
FTG limit due to the fermionic statistics. The half-width 
of the profiles become larger as the attractive interaction 
increases. 

As a comparison, the momentum distributions of 
Bosons with repulsive interactions are given in Fig. 4 
which are obtained by the Fourier transformation of the 
one-body density matrix shown in Fig. 5. Obviously, 
the momentum distribution for the bosons exhibits quite 
different behaviors. As shown in Fig. 4, there is an ob- 
vious peak around the zero momentum point and the 
height of the peak shrinks with the increase of the repul- 
sive interaction. Furthermore we find no oscillation at all 
in the momentum distribution for bosons. Even in the 
Tonks limit, the momentum distribution of Bosons does 
not show shell structure like free Fermions. The largest 
probability of distribution appears around the zero point 
of momentum and decreases rapidly for finite momen- 
tum values. Stronger repulsive interaction between the 
bosonic atoms tends to spread out the distribution into 



higher momentum space. Although the momentum dis- 
tribution for the interacting bosonnic gas es has been 
studied by different numerical methods [35l |3(| , an ex- 
act result has never been given except in the TG limit 

From the eqs. ([7]) and ([8]), we see that the energy 
level structure of our fermionic model is exactly the 
same as the corresponding bosonic model with inter- 
action strength related by c = — 1/cp. Therefore the 
thermodynamic properties of the fermionic atoms with 
attractive p-wave interactions are the same as the well 
known thermodynamic properties of the ID boson gas 
[29l l33l . l34j with inverse coupling c = —1/cp. This im- 
plies that there exists a Bose-Fermi duality between the 
p-wave fermionic model and bosonic one with arbitrary 
interactions which can not be distinguished by the ther- 
modynamic properties. However, due to the different ex- 
change symmetry of the wave functions, the observables 
associated with the wave functions rather than the square 
of wave functions (density distribution) should display 
different behaviors. As we have shown, the off-diagonal 
density matrix and the momentum distributions are dif- 
ferent greatly, which result from the different statistics 
properties followed by Bosons and Fermions. 



IV. CONCLUSIONS 

In summary, we have investigated the ground-state 
properties of fermionic gases with attractive p-wave in- 
teractions in a one-dimensional hard-wall trap. With the 
Bethe ansatz method, the explicit wave function of the 
ground state and therefore the one-body density and mo- 
mentum distributions are obtained. It turns out that the 
density distributions show one-to-one correspondence be- 
tween Fermions with attractive p-wave interactions and 
Bosons with repulsive interactions. For weak attrac- 
tive interaction the density distributions of Fermions dis- 
play shell structures and the Boson-like distributions ap- 
pear as the interactions become stronger. In the FTG 
limit, the Fermi gas exhibits the same distribution as 
that of the ideal Bose gas. This again confirms the 
Bose-Fermi duality: strongly interacting Bosons behave 
like Fermions, and vice versa. Nevertheless, from the 
viewpoint of momentum distribution, the conclusion is 
rather different. In the full interacting regime the mo- 
mentum distributions of Fermions show typical oscilla- 
tions, which is in sharp contrast with the case of bosonic 
atoms. In the Tonks limit of infinite interaction, although 
the density profiles of one dimensional bosons display the 
Fermion-like distribution, the momentum distribution is 
still Boson-like. 
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